Channel plasmon-polaritons: modal shape, dispersion, and losses 
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We theoretically study channel plasmon-polaritons (CPPs) with a geometry similar to that in recent 
experiments at telecom wavelengths (Bozhevolnyi et al, Nature 440, 508 (2006)). The CPP modal shape, 
dispersion relation, and losses are simulated using the multiple multipole method and the finite difference 
time domain technique. It is shown that, with the increase of the wavelength, the fundamental CPP mode 
shifts progressively towards the groove opening, ceasing to be guided at the groove bottom and becoming 
hybridized with wedge plasmon-polaritons running along the groove edges, (c) 2008 Optical Society of America 
OCIS codes: 240.6680, 130.2790, 260.3910. 



The guiding of light within a subwavelength cross sec- 
tion has been recently attracting a great deal of atten- 
tion due to ever increasing demands for miniaturization 
of photonic circuits. Light may be confined in the di- 
rection perpendicular to a flat metallic surface for ener- 
gies below the metal plasma frequency The mode guided 
along the metallic interface is known as surface plasmon- 
polariton (SPP). Various geometries have been proposed 
to achieve confinement of the plasmon-polariton in the 
plane transverse to the propagation direction. Among 
these proposals, the plasmon-polariton guided by a V- 
shaped groove carved in a metal (channel plasmon- 
polariton, CPP) is particularly interesting. CPPs were 
theoretically suggested by Maradudin and coworkers^ 
and subsequently studied in the visible regime.^' ^ Re- 
cently, CPPs have been experimentally investigated at 
telecom wavelengths,^ displaying strong confinement, 
low damping, and robustness against channel bending. 
Thank to these properties, prototypes of basic devices 
could be demonstrated.^ The mentioned devices have 
been developed with the help of the effective index ap- 
proximation but, to our knowledge, no rigorous electro- 
dynamic computation of CPPs at telecom wavelengths 
has been reported. The effective index approximation 
can deliver information about the dispersion relation, 
but it is expected to be inaccurate for frequencies close 
to the mode cutoff and is unable to determine modal 
shape and polarization. The functionality of many de- 
vices relies on the overlapping of electromagnetic fields at 
various sites inside the device. For this reason the knowl- 
edge of the modal shape is essential to provide a solid 
foundation for the design of CPP-based devices. Here 
we present rigorous simulations of guided CPPs aimed 
to elucidate their characteristics at telecom wavelengths, 
including full vectorial modes, dispersion, and losses. We 



show that, contrary to what is commonly believed, CPPs 
at telecom wavelengths are not guided at the groove bot- 
tom, at least for the groove parameters used in the exper- 
iments.^'^ Instead, the CPP field at the groove entrance 
hybridizes with wedge plasmon-polaritons (WPPs) run- 
ning along the edges of the groove. 

Our goal is to understand the fundamental CPP mode 
guided by realistic grooves at telecom wavelengths.^ Nev- 
ertheless, in order to comprehend the behavior in this 
regime, which is close to cutoff, we will consider a broader 
spectrum, higher order modes, and a number of different 
geometries. The simulations have been performed with 
two rigorous electrodynamic techniques: the multiple 
multipole method (MMP)^^ and, where mentioned, the 
finite difference time domain (FDTD) method. Within 
the MMP method the corners are rounded (10 nm radius 
of curvature). FDTD results were converged for a mesh 
of about 5 nm. Such fine meshes are essential, the more 
so for wavelengths shorter than 0.8 /am. The grooves are 
carved in gold and we employ experimentally measured 
values of the dielectric permittivity e. 

Figure 1 shows the dispersion relation for a non- 
truncated groove with a semiangle of = 12.5° and in- 
finitely long sides. This structure sustains two modes, 
being termed CPP(oo) (see right insets), which are out- 
side the dispersion line of the SPP at a flat surface. The 
modal shape (time averaged electric field) is shown in 
the right insets for a wavelength of A = 0.6 /j,m. In the 
figure it is also plotted the dispersion relation for a non- 
truncated metallic wedge of semiangle a = 51.25° and 
infinitely long sides. The corresponding wedge mode run- 
ning along the edge is termed WPP(oo) (see left inset). 
WPP(oo) for this a will be relevant when we later trun- 
cate the above groove at a finite height: it corresponds to 
the edges at both sides of the finite- height groove. The 
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Fig. 1. (Color online) Dispersion relation for various 
modes. Black thick line: SPP mode on a flat surface. 
Green lines (squares): CPP(oo) modes for an infinitely 
deep groove. Red line (open circles): WPP(oo) mode for 
an infinitely deep wedge. Right insets: time averaged 
electric field of the two CPP(oo) modes at 0.6 /im. Left 
inset: same for the WPP(oo) mode. The lateral size of 
the insets is 2 //m. 

WPP(oo) modal field at 0.6 /im is shown in the left inset. 
For increasing wavelength all three modes approach the 
SPP line (none of them has a cutoff). In this process the 
modal shapes remain qualitatively the same, the only 
difference being that the fields are expelled away from 
the groove or wedge corners. 

Figure 2 represents a similar plot but now a groove 
of finite height is considered, the height being 1.172 /im. 
The CPP modes exhibit now cutoff at different wave- 
lengths (~ 1.44 /im for the first mode and ~ 0.82 /im 
for the second one). This was advanced in Ref. 12, and 
it is a consequence of the above mentioned behavior of 
the fields for increasing wavelength. As the wavelength 
grows, the field is pushed out of the groove and, after 
a certain threshold, it can no longer be confined by the 
groove sides and is radiated in the form of SPPs along 
the contiguous horizontal metal surfaces. It is important 
to realize that, before reaching the SPP dispersion line, 
both modes approach and cross the WPP(oo) line. This 
means that close to cutoff the CPP modes must be hy- 
bridized with the modes running on the edges at both 
sides of the groove. This idea is visualized in the in- 
sets, that render the modal shapes (time averaged elec- 
tric field) at 0.6 /im. At this wavelength the first mode 
is not close to WPP(oo) and the hybridization does not 
take place, but it is already happening for the second 
mode. The described phenomenon is even more distinct 
in Fig. 3 displaying the fundamental mode for increasing 
wavelengths. It is observed that the CPP mode becomes 
more and more mixed with the WPP(oo). Close to cut- 
off (at about 1.44 /im) the mode is not guided at the 
groove bottom anymore but rather at the groove edges. 
A hint of this possibility was mentioned in Ref. 13. In 
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Fig. 2. (Color online) Dispersion relation for various 
modes. Black thick line: SPP mode on a flat surface. 
Blue lines (full circles): CPP modes for a groove of 
height 1.172 /im (computed with MMP method). Trian- 
gles: same as before computed with FDTD method. Red 
line (open circles): WPP(oo) mode for an infinitely deep 
wedge. Insets: time averaged electric field of the two CPP 
modes at 0.6 /im. The lateral size of the insets is 2 /im. 

the experiments, the edges at both sides of the groove 
have larger radius of curvature than in the previously 
presented simulations. We have verified that this does 
not alter our conclusion by repeating the same compu- 
tation with a radius of curvature of 100 nm at the groove 
edges (and keeping 10 nm at the bottom). Figure 3(d) 
shows the instantaneous transverse electric field for this 
case and it is clear that hybridization with edge modes 
still occurs. The transverse electric field is approximately 
horizontal inside the channel (an assumption used by 
the effective index approximation), but it is not horizon- 
tal near the edges where the field is maximum. Let us 
note in passing the excellent agreement of the two tech- 
niques employed here (the residual discrepancy in Fig. 2 
for the fundamental mode at 0.6 /im is due to different 
rounding schemes of the groove bottom in the two meth- 
ods). From the point of view of fabrication it is useful 
to mention that, for A G (0.6 /im, 0.8 /im), the disper- 
sion relation is extremely sensitive to the fine details of 
the groove bottom (e.g., rounding), as concluded after a 
large number of simulations where the details of the bot- 
tom were subjected to small perturbations. On the other 
hand, this does not happen for telecom wavelengths (as 
expected from the modal shape), a circumstance that 
has also been observed experimentally.^ Note that the 
calculated cutoff wavelength of the fundamental mode is 
somewhat lower than the wavelengths used in the experi- 
ments. This discrepancy can be ascribed to (small) differ- 
ences in the groove geometry, both in the groove shape 
(angle, side flatness) and in the groove depth, and/or 
different dielectric permittivity of gold. We have verified 
(not shown here for brevity) that slightly less negative £ 
or / and smaller groove semiangle 9 leads to a higher cut- 
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Fig. 3. (Color online) Modal shape of the CPP fun- 
damental mode for increasing wavelength A. (a) A = 
0.6 /im, (b) A = 1 /im, (c) A = 1.4 /im (close to cutoff). 
These panels display the time averaged electric field, (d) 
Instantaneous transverse electric field at A = 1.4 /im for 
a structure with groove edges rounded with 100 nm ra- 
dius of curvature. All panels have a lateral size of 2 /xm. 



off wavelength. Finally, the experiments were conducted 
at ambient conditions so that water condensation could 
not be excluded (a very thin water layer can significantly 
increase the cutoff wavelength). 

The effect of absorption is summarized in Fig. 4 that 
renders the propagation length / = [2lm{kz)]~^ , versus 
wavelength for the various considered structures {kz is 
the modal wave vector). The propagation lengths are in 
all cases smaller than that of SPPs at a fiat surface. This 
is a consequence of the field enhancement at the cor- 
ners and the field confinement that decreases the por- 
tion of field propagating in air. When comparing the 
CPP modes it is observed that the effect of truncation 
at a finite height is only important for wavelengths larger 
than 1 /im, which is reasonable because the field is very 
much confined at the groove bottom for smaller A. For 
longer wavelengths the CPP propagation length is de- 
creased as compared to that of CPP(oo). At A = 1.4 /im 
we find that /cpp = 53 /im. The values reported in Ref. 8 
at 1.55 /im are twice as large. The discrepancy can be 
again ascribed to slight differences in geometry and/or 
dielectric permittivity that rise the cutoff wavelength. If 
the trend of the line corresponding to the CPP is ex- 
trapolated, we find good agreement with the reported 
data. It is to be observed that the propagation length of 
WPP((X)) is significantly higher, a fact that could find 
obvious applications. 

In conclusion, we have presented rigorous computer 
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Fig. 4. (Color online) Propagation length versus wave- 
length for various modes. Black thick line: SPP mode 
on a fiat surface. Blue line (full circles): CPP fundamen- 
tal mode for a groove of height 1.172 /im. Green line 
(squares): CPP((X)) fundamental mode for an infinitely 
deep groove. Red line (open circles): WPP(oo) mode for 
an infinitely deep wedge. 

simulations of CPPs at telecom wavelength. CPPs have 
been fully characterized in terms of modal shape, dis- 
persion, and losses. We have shown that, for relatively 
shallow grooves, the field is guided at the groove opening 
and is hybridized with modes running along the groove 
edges (WPPs). We expect that our findings will be of 
help for the design of improved CPP devices. 
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